Single particle spectrum of resonant population imbalanced Fermi gases 
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We use a T-matrix approximation to calculate the single particle spectrum of the normal state of 
a gas of Fermionic atoms at low temperature. In the strongly interacting regime of the polarized gas, 
we find that the spectrum is separated in two branches, leading to a double-peaked radiofrequency 
spectral feature. 
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The experimental observation of fermionic superfluid- 
ity in ultracold alkali gases [if is one of the most impor- 
tant achievements in atomicphysics in recent years. By 
applying a magnetic field [3, [Sj, experimentalists have 
been able to control the interaction strength and ob- 
served how superfluidity evolves as one increases the at- 
traction between particles from the weak-coupling BCS 
limit. As was predicted by theory, they saw that the su- 
perfluid state of weakly bound pairs continuously evolves 
to one described as a Bose Einstein condensate (BEC) of 
tightly bound molecules Q- In the intermediate region 
the scattering amplitude saturates the bounds set by the 
unitarity of the S-matrix, and a strongly interacting su- 
perfluid with universal properties appears. The mecha- 
nism at the root of fermionic superfluidity is the pairing 
between atoms in two different states, which we will refer 
to as t and J, spin. In recent experiments 0, where 
the spin relaxation time exceeds all other experimen- 
tal timescales, researchers have created spin imbalanced 
Fermi gases, where the ratio of the number of atoms in 
each spin state Ni/N^ = a; is below one. This popula- 
tion imbalance suppresses superfluidity 0, 0], allowing 
one to study a low temperature normal phase. Unlike 
the superfluid phase, where symmetries determine most 
of the systems behavior, the normal phase could be quite 
exotic. Here we use a T-matrix approximation to study 
the single particle spectrum of the normal state of an 
imbalanced Fermi gas, calculating thermodynamic prop- 
erties and experimental observables. 

We will be particularly focused on radio frequency 
(RF) spectroscopy, where one uses radio waves to drive 
a transition from one of the two spin states (say J,) to a 
third, excited state |ex). By exploring the rate of trans- 
fer as a function of the frequency of the radio waves, 
one can gain information about the many-body density 
matrix. This technique was used to measure the bind- 
ing energy of pairs in the BEC regime P , and extended 
across the resonance to learn about the pairing gap of uni- 
tary gas [l^ . Recent measurements on imbalanced gases 
of ^Li atoms have shown the remarkable feature that the 
RF spectroscopy is qualitatively unchanged by polarizing 
the gas Uj. At high temperature (in the normal state 
of either the equal spin gas or the spin imbalanced gas) 
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FIG. 1: Three dimensional depiction of mean field spectral 
density of an unpolarized superfiuid gas: (a) "BCS regime", 
A = 0.5/i, (b) "BEC regime" A = -10/x. The points are 
raised a distance above the k-u plane proportional to the 
height of the delta-function peaks. These large values of A 
are chosen to emphasize the important qualitative features. 



one sees a single peak in the RF spectrum. As the tem- 
perature is lowered a second peak appears. In both the 
unpolarized and spin imbalanced cases the appearance 
of this second peak occurs at a temperature of order the 
superfluid transition temperature of the gas with equal 
spin populations. At lower temperatures yet, the original 
peak disappears. These results are surprising. Although 
the unpolarized gas undergoes a superfluid phase tran- 
sition, the spin imbalanced gas remains in the normal 
state. How can the normal state of the imbalanced gas 
have an RF spectrum which is qualitatively indistinguish- 
able from that of the superfluid state of the equal spin 
gas? Is the ground state of the spin imbalanced gas an 
exotic state containing noncondensed bosonic pairs? 

As we show below, the RF spectra observed by Schunck 
et al. [lH are consistent with a traditional Fermi liquid 
scenario, where the low energy excitations can be placed 
in one-to-one correspondence with those of an ideal Fermi 
gas. The higher energy excitations of this system are 
however quite unusual. We find an anomaly in the down- 
spin spectral density at energies/momenta around the 
up-spin Fermi energy/momentum. This spectral feature, 
characterized by a sudden bending of the spectrum, ef- 
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fectively separates the single particle spectruni into two 
branches, yielding the two peaks seen in experiments. 

Before discussing a theory of the spin imbalanced nor- 
mal state it is helpful to review the physics of RF spec- 
troscopy in mean-field theory of the superfluid state. For 
simplicity we will neglect the interaction between the 
atoms in the |ex) state and the other two spin states. 
This approximation should not be quantitatively correct 
for the recent experiments, where there are several Fesh- 
bach resonances in close proximity Despite the im- 
portance of understanding these final state effects 
one does not expect that they play any qualitative role 
in the observations. Under these assumptions, parti- 
cles in the excited internal state have a free dispersion 
eex(fc) = eex(O) -I- fc^/2m, and Fermi's golden rule tells 
us that the number of atoms transferee ! by an RF spec- 
troscopy experiment is proportional to |14l . i 
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where v is the detuning of the radio frequency field from 
the free-space splitting between the |ex) state and the [ 
state, ri(e) = (e^^ -I- 1)^^ is the Fermi function, /3 is the 
inverse temperature, and the spectral density Ai{k,u)) 
represents how many down-spin single particle states ex- 
ist at a given momentum and energy (measured from the 
down-spin chemical potential). Figure [T] shows the mean- 
field result Aik^uj) — vl6{u! — Ek), where the coherence 
factor is vl = [e^ — Ek)/{2Ek), the noninteracting spec- 
trum is efe = fc^/ (2m) — /x and the quasiparticle spectrum 
is El = el + A"^, with the superfluid gap given by A. 
The spectral density is qualitatively different in the BCS 
regime (/i > 0) and the EEC one < 0). Qualita- 
tively, the unitary gas is similar to the BCS limit as it 
has /i > 0. It is important to note that aside from near 
the places where the two branches become close, most of 
the spectral weight lies near the free particle spectrum 
uj = efc. 

These mean-field spectra can be understood by noting 
that within the mean-field picture there are two ways to 
add an down-spin particle with momentum k. One can 
directly add it (costing energy e^i) or one can remove 
an up-spin particle of momentum q — k and simultane- 
ously add a pair with momentum q. In this latter case 
the change in energy is the difference between the energy 
of the pair measured from the pair chemical potential 
Eb + q^/2mb — fib, and the energy of the up-down parti- 
cle which was removed. In the superfluid state there is 
a condensate of g = pairs, requiring E^ — fib = 0. One 
therefore finds that the energy to add the spin-down par- 
ticle in this manner costs energy — e^f . An avoided cross- 
ing between these two modes gives the BCS spectrum. In 
the BEC limit there is no avoided crossing, and one just 
has two parabolic bands. 

In contrast to the mean field treatment of the super- 
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FIG. 2: Spectral density of spin down quasiparticles A{e, k) 
at zero temperature. The e axis is normalized by pi, the k 
axis by ^J2m[s,^ , m being the bare mass of atoms. The value 
of the chemical potential ratio in this plot is /xj.//X| = —0.6, 
the scattering length is taken to be infinite. 



fluid, when one includes fluctuations in the normal state 
(either above Tc or in the normal state of a spin imbal- 
anced gas), there are many possible channels for adding 
a single down-spin particle. In principle, one should con- 
sider all possible ways of exciting the many body state 
while adding the particle. Below we will take into account 
processes where a single spin-up particle hole pair is cre- 
ated. Including all such processes leads to anomalous 
features in the spectral density of spin-down particles. 
The spectral density is a nonanalytic function of energy 
and momenta at the values where the phase space for the 
creation of particle-hole pairs hits the boundary of the 
particle-hole continuum. This feature, which is a conse- 
quence of the sharp Fermi Surface, corresponds to Kohn 
anomalies which are observed in the spectra of phonons 
[r?! . At momenta of the order of the Fermi momentum 
of spin up particles, the spectrum of spin down particles 
bends and becomes damped, creating a dip in the density 
of states. This effectively translates into a separation of 
the spectrum into a low momentum (shifted down from 
the free spectrum) and a high momentum branch (shifted 
up from the free spectrum) , which gives to the density of 
states in the normal state (Fig. 2) a structure similar to 
that in the superfluid state (Fig. 1). 

Geometrically, the RF spectroscopy intensity can 
be found by overlaying the parabola k'^/2m—v-\-fii on the 
down-spin spectral density graph and quantifying their 
overlap. Note that since momentum is conserved under 
absorption of an RF photon, this experiment is very dif- 
ferent from a point-contact tunneling experiment in solid 
state physics. Rather it is equivalent to a solid state 
photoemission experiment (if one does not momentum- 
resolve the final state) or a tunneling experiment in which 
momentum is conserved (such as occurs in tunneling be- 
tween parallel wires [18]). 



3 




FIG. 3: Single particle density of states P|(e), normalized by 
fj,1 , plotted in the absence (dashed) and in the presence (solid) 
of interaction as a function of e/pi|. The density of states is 
reduced relative to the noninteracting value near e ~ /i| . 



In the superfluid state at finite temperature one sees 
a bimodal RF spectrum as each of the two branches of 
the single particle spectrum contributes to the RF hne- 
shape at different frequencies. Since the upper branch 
asymptotically approaches that of a free gas, one of the 
two spectral lines begins at v = 0. The separation S be- 
tween the two lines is given by the difference between the 
k — energy of the free parabola and the fc = energy 
of the lower branch: 6 - 



+ — (It is particu- 
larly important to realize that this splitting is not simply 
the gap A [3].) At high temperature the contribution 
from the upper branch dominates, while at low temper- 
ature the contribution from the lower branch dominates. 
On the BCS side of resonance, where pairing mostly oc- 
curs near the Fermi surface, both peaks are dominated 
by unpaired states. As is clear from the spectral density 
shown in Figure [ija), the presence of pairs only plays a 
role very near to the gap. On the other hand, on the 
BEG side of resonance, the entire lower branch is due to 
pairs. Clearly, the existence of a shifted low temperature 
RF peak is not sufficient to conclude that the system is 
superfluid, a counterexample being a gas of fermions with 
interaction treated in the Hartree Fock approximation. 

To extend this picture to the spin imbalanced gas we 
must calculate the spectral density A{k,uj). Since there 
exists to our knowledge no controlled method for calcu- 
lating the properties of a Fermi gas at unitarity, we work 
with the simplest approximation which captures the basic 
physics - namely the many-body T-matrix approxima- 
tion which generalizes Noziere and Schmidt- Rinks theory 
of the BCS-BEC crossover to the spin imbalanced gas. 
Variants of this approximation have been used by several 
groups 0,j2l|. We use the one described by Combescot 
et al 22, l23|, who have shown that this approxima- 
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FIG. 4: Calculated RF spectrum in the normal state at uni- 
tarity: (a) T = (b) r = 2^1, (c) T = (d) T = 0. Solid 
lines correspond to [iij [i\ = —0.60671, resulting in vanishing 
minority species density at T = 0, Dashed lines correspond 
to /ii//iT = 0.1. 



tion gives remarkably good agreement with Monte-Carlo 
calculations in the limit of vanishing downspin density. 
Specifically we take 



^X(tj, k) = Im E4(tJ, fc)j 

I]x(c^,fc)= / dzrx(z,A:)/(24c^-^)) 



(2) 
(3) 



rx(c.,fc)= / A/(27r)3A(c^-f ej,fc-f g) (4) 

A(w, fc) = 2Imr(a;, fc) (5) 
r(tj,fc) = (47r;iVm)/(a"^ + (6) 



] 27r(w - z) j {2tiY 
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where e% = k'^ /2m - fi^, and = 0{-el). Note that 
spectra used to calculate the self-energy are free spec- 
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tra. Self-consistency should not qualitatively change the 
spectral density, and may actually make the theory less 
accurate 
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The resulting spectral density is shown in Fig. [21 where 
the separation between a small momentum branch, un- 
derdamped and parabolic, and a branch at higher mo- 
menta, shifted upward, is clearly recognizable. This two 
branch structure is similar to the mean-field spectral den- 
sities in Fig. [U except for the fact that the branches are 
joined by an overdamped continuum. In the limit of van- 
ishing Ni, (occuring at = —0.6067) we find that the 
lower branch is described by E{k) — fc^/2m*, with an 
effective mass m* = 1.16. At sufficiently large momenta, 
the damping of quasiparticles is vanishingly small. This 
is because the fundamental interaction is short-ranged, 
and unitarity of the S-matrix requires that scattering 
becomes weak at large momenta. Figure [3] shows the 
integrated density of states 
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where a dip is clearly visible at energies close to the Fermi 
level. Here, we mention that a similar dip is found in the 
spectral density of cuprate superconductors ^20J . 

In Figure 4, we use Eq. 1 to calculate the RF- 
spectroscopy lineshape. Note that since we are using the 
zero-temperature spectral density the finite temperature 
line-shapes are at most qualitative. The general struc- 
ture, however is generic. At the lowest temperatures only 
the bottom branch of the spectrum is occupied, and one 
sees only a single peak, shifted from (5 = by an amount 
proportional to /i^. In the limit of vanishing downspin 
density this shift is directly equal to the downspin chemi- 
cal potential /i| — = 0, = 0), and provides a model 
independent way to determine this quantity. As temper- 
ature rises the upper branch becomes occupied, resulting 
in a second peak at lower frequencies. The weight in 
the second peak grows with the temperature, eventually 
dwarfing the low temperature peak. 

We believe that in a finite temperature calculation of 
A{k,u!) one would find that the separation between the 
two branches would become smaller as temperature in- 
creased. This would manifest itself in the low tempera- 
ture peak slowly moving towards lower energy, merging 
with the high temperature peak. The splitting should 
vanish at a characteristic temperature T* ~ pL^/kg. 
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